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Abstract 

A two-dimensional square-lattice model for the formation of secondary structures in proteins, the 
hydrogen-bonding model, is extended to include the effects of solvent quality. This is achieved by 
allowing configuration-dependent nearest-neighbour interactions. The phase diagram is presented, 
and found to have a much richer variety of phases than either the pure hydrogen-bonding self- 
avoiding walk model or the standard ©-point model. 
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I. INTRODUCTION 



Self-avoiding walk models have been used for many years as models of real polymers in 



solution 171 ] . The thermodynamical behaviour of a linear polymer in a good dilute 
solvent is dominated by its entropy, which may be well modelled by an excluded volume 
interaction, leading to the idealised model of a self-avoiding walk on a lattice. As the tem- 
perature of the polymer is lowered, typically the quality of the solvent is degraded, and 
the difference in affinity between the monomers (chemical building blocks of the polymer) 
and between the monomers and the solvent molecules becomes important. At low enough 
temperatures the polymer collapses and will precipitate from solution. This difference in 
affinities may be modelled in the self-avoiding walk model by an effective attractive interac- 
tion between neighbouring steps of the walk. The high temperature (good solvent) and low 
temperature (bad solvent) regimes are separated by a phase transition point known as the 
6-point. 

The canonical model for this system is the O-point model, which consists of placing the in- 
teractions between nearest- neig hbour lattice sites which have been visited non-consecutively 
by the self-avoiding walk[?], Typically this model is studied in the grand-canonical en- 
semble where the length of the walk is governed by a step fugacity K. The grand-canonical 
partition function, Z, is then written: 

Z= K n t n ' } (1) 

walks 

where N is the length of the walk, Nj are the number of nearest-neighbour interactions and 
r = exp(—j3e), j3 = 1/kT, e < is the (attractive) energy gained per nearest-neighbour 
contact. This model describes well the behaviour of simple linear polymers in solution. The 
phase diagram in the fugacity/temperature plane is shown in figure [TJ It is expected that 
any typical size of the walk, such as the radius of gyration, should scale as a power law with 
the length of the walk as follows: 

Rg ~ (N) v , (2) 

where v is simply the correlation length exponent defined in magnetic models. Indeed, the 
self- avoiding walk model may be mapped onto an 0(n) symmetric spin model in the limit 
n — > oj^j]. It may then be seen that the phase transition line shown in figure [1] is the line 
on which the average length diverges, for large temperatures (small (3) continuously and for 




FIG. 1: The phase diagram for the standard 0-point model, showing a low- if zero density (finite 
walk length) phase and a high-ET (critical) collapsed phase, where the walk density is finite. At 
high temperatures (low (3) the transition is second order (solid line), whilst at low temperatures 
(high P) the transition becomes first order (dashed line). These two behaviours are separated by 
a tricritical point, the O point. 

low temperatures discontinuously. The O-point is then identified with a tri-critical point. 

Whilst this model is of most practical importance in three dimensions, it has been exten- 
sively studied in two dimensions, because it is expected that the critical behaviour is richer 
in two dimensions, particularly since the upper critical dimension of a tri-critical point is 3. 

An interesting question arises: the basic ingredients modelled by the G-point model 
which enable it to capture the essence of the real polymer in solution are the modelling of 
the entropic repulsion by the excluded volume interaction (self-avoidance) and the modelling 
of the short ranged attraction between the monomers. It should be expected, then, that any 
consistent way of modelling these two features would lead to a thermodynamically equivalent 
model. With this in mind, Blote and Nienhuis J] introduced a variant on the model in which 
the self-avoidance restriction is relaxed in that the walk can now visit lattice sites more than 
once, however the walk is not allowed to visit the lattice bonds more than once, and the 
walk is not allowed to cross itself. The interactions are now introduced for the doubly visited 
sites. Surprisingly this model has a phase diagram which differs from the O-point model in 
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many important respects: there is an extra phase transition in the dense region of the phase 
diagram. This phase line is in the Ising universality class with a value of the correlation 
length exponent v = l^, [ijj]. The new collapse transition is not in the same universality 
class as the point, having an exponent v = 12/230, as compared to v — 4/7. At first 
sight these differences seem to be in contradiction with the idea of universality, which is 
required if we are even to think of modelling a polymer in solution by a lattice based walk 
model. In fact, Universality is not violated. If the walk fills the lattice with a finite density, 
and its fractal (Haussdorf) dimension is the same as the lattice dimension, then the walk 
sees the lattice, and may be subject to lattice frustration effects. At the point the density 
is zero, whilst it was shown that the density at the collapse transition in the Nienhuis-Blote 
model is non-zero[lo|. The presence of an Ising transition in the dense region is an indication 
that the lattice interactions tend to pick out a preferred direction, here corresponding to the 
lattice diagonals. 

Other lattice models have been introduced which contain collapse transitions, notably 

n n n 

the bond- interacting self-avoiding walk[3|, 18|, H6( and the Hydrogen-bonding self- avoiding 
walk[l|, The first is simply the G point model in which the interactions are now between 
the nearest-neighbour visited lattice bonds. The Hydrogen-bonding self-avoiding walk was 
introduced to model the formation of secondary structures in proteins under the influence 
of the hydrogen-bond. Hydrogen bonds are induced by dipole-dipole interactions, and im- 
pose an orientation on the interacting portions of the polymer. How these interactions are 
implemented in the Hydrogen model is shown in figure EJ 

The bond-interacting self-avoiding walk has been studied using mean-field type calcula- 
tions on the Bethe lattice [3J] and the Husimi lattice [l6|. These different studies have proposed 
radically different phase diagrams. A recent transfer matrix calculation 8| indicates that the 
correct phase diagram is that proposed by Buzano and Pretti|3], and shown schematically in 
figure [3l Whilst there is a collapse transition in the same universality class as the standard 
B point, there is also the presence of an additional transition line in the dense region of the 
phase diagram. Unlike the Nienhuis-Blote model, this transition line seems to be exotic in 
nature, with a non-divergent susceptibility Q]. 

The Hydrogen-bonding model was also investigated using transfer matrices 111]. The 
phase diagram was schematically similar to the Nienhuis-Blote model, but the collapse 
transition was found to be first order. The high density transition line seemed not, in this 
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(a) 



(b) 



(c) (d) 

FIG. 2: The nearest-neighbour interactions are split into two classes, those of type (a) where 
four bonds forming two parallel lines model the hydrogen bonds, whilst the others (b,c and d) 
model the solvent interactions. Configuration (a) induces a preferred orientation, whilst the other 
configurations do not. 

study, to be in the Ising universality class, but the error bars were so large as to make it 
hard to rule out any possibility. We return to this question in this work, and using recently 
developed corner-transfer-matrix renormalisation-group (CTMRG) methods[9] we manage 
to give an accurate estimate of the critical exponent u, clearly ruling out any possibility that 
it could be in the Ising universality class. 

Faced with this variety of different behaviours, it almost seems that any change to the 
model leads to different behaviour for dense interacting self-avoiding walks. In order to 
investigate the relationship between these different behaviours we introduce an extension of 
the Hydrogen-bonding self-avoiding walk model to include 0-type interactions. These inter- 
actions are shown in figure [2J There is also a more direct motivation for these interactions; 
whilst a protein folds under the influence of the hydrogen bonds, it is still subject to the 
quality of the solvent. Indeed it is the quality of the solvent (or physiological temperature) 
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FIG. 3: A schematic version of the phase diagram for the bond- interacting self-avoiding walk 
model, proposed by Buzano and Pretti[3j and confirmed by Foster [gl . Phase I is the low-K finite 
walk phase, II is the critical collapsed (liquid) phase and III is the crystalline oriented phase. 

which decides whether or not a protein is functional. We will show that this enlarged model 
contains several, if not all, of the different thermodynamic behaviours found above, and 
provides a unifying framework in which to understand the effect of geometrical frustration 
in lattice walk models. 

In the next section the model is presented. In section II II I results are first presented for 
two cases, showing two different behaviours. The results found then enable a mapping of 
the full phase diagram. The article ends with a discussion of the results. 

II. THE MODEL 

The model studied in this article involves the self-avoiding walk on the square lattice with 
interactions between non-consecutive visited nearest-neighbour sites on the lattice. Unlike 
the standard G-point model, the interactions are split into two sets, as shown in figure [2, 
between those which specify a particular direction, the hydrogen bonds, and those that do 
not, the solvent interactions. Hydrogen bonds carry an interaction energy —eh and the 
others carry an interaction energy — e. The thermodynamic behaviour may be investigated 
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by introducing the grand-canonical partition function, 2, from which many of the relevant 
thermodynamic quantities may be calculated. The grand-canonical partition function is 
given by: 

Z = K N exp (p (Nje + AW)) , (3) 

walks 

where Ni are the number of solvent interactions, and Nh are the number of hydrogen bonds. 
The fugacity, which controls the average length of the walk, is denoted by K, and A" is the 
total length of the walk. For convenience we define a = e/eh, and without changing the 
physics of the model, we may set £# = 1; this simply sets the temperature scale. The 
partition function then becomes: 

Z=Y. RN ex P iP ( N h + Nra)) . (4) 

walks 

The primary tool we use in this article to map out the phase diagram of this model is the 
transfer matrix. The transfer matrix method involves studying the model on a lattice strip 
of length L x — > oo and width L y , finite. In its simplest form, the model is considered on a 
lattice with periodic boundary conditions in both the x and y lattice directions. In this case 
the partition function may be rewritten in terms of a matrix product: 

Z = TyT l *, (5) 

where T is the transfer matrix, which contains all the factors required to take account of all 
possible walk configurations between any two given lattice columns. Details on the transfer 
matrix method used in this article may be found in reference la. 

The partition function may then be expressed in terms of the eigenvalues A», of the 
transfer matrix T: 

z = J2 x f x - ( 6 ) 

i 

The dimensionless free energy per spin is given by 

f = -^— log Z. (7) 

In general the largest eigenvalue is non-degenerate, and the sum is dominated by this 
largest eigenvalue, Ao, giving, in the limit L x — > oo, 

/ = J_l ogAo . ( 8 ) 

Ly 



The problem is now reduced to studying the behaviour of the thermodynamic quantities 
as a function of the width, notably using finite-size scaling techniques. In what follows we 
will drop the subscript y and denote the lattice width by L. Once the free energy has been 
calculated, other quantities of interest can be calculated by taking suitable derivatives, for 
example the density of the walk on the lattice is given by 

It is however possible to calculate such quantities directly from the eigenvalues and eigen- 
vectors of the transfer matrix. To see this, it is necessary first to calculate the probability 
of having a given walk configuration C x in column x. This probability is simply the ratio of 
the partition function restricted to having configuration C x in column x and the unrestricted 
partition function, which in terms of transfer matrices may be written: 



Tr {T*\C,){C,\T L -- X } 



*» = i^lfT5 (10) 

Writing \C X ) in terms of the eigenvectors, \i) of T gives: 

m = lim ^xtmm 

p(C) = (0|C) 2 , (12) 

where the eigenvectors are normalised. The subscript x may be omitted by invoking trans- 
lation invariance. The density, for example, is then found using 

c L 

= E^r 1 (oi0 2 , (is) 

c L 

where N(C) is the number of occupied lattice bonds in configuration C. The susceptibility 
can then be calculated either by taking a derivative of the density, or by calculating directly 
(N 2 ) for the column, and hence the fluctuation. The two methods give slightly different 
results for a finite width strip, but agree in the thermodynamic limit. In the present article 
we choose to calculate the fluctuation directly. 

It is straightforward to show that the correlation length, £, depends on the largest two 
eigenvalues through 



If the two eigenvalues become equal in modulus, the correlation length diverges, which is 
characteristic of the long-range order found at a critical point. For an integer spin model 
(Ising, XY, Heisenberg etc.), the transfer matrix is positive (all elements strictly larger than 
zero) and Frobenius' theorem states that the largest eigenvalue is non-degenerate for finite 
matrices. This implies that the correlation length may only diverge in the thermodynamic 
limit L x , L y — > oo. In our case, however, the transfer matrix is sparse, and may be block 
diagonalised into an odd and an even sub-matrix. The odd sub-matrix is the transfer matrix 
for the walks which cross the lattice in the x-direction an odd number of times, whilst the 
even sub-matrix is the transfer matrix of walks which cross the lattice an even number of 
times. We include in the even sub-matrix the empty lattice configuration. There is no 
mathematical reason why the largest eigenvalues of the different sub-matrices should not 
coincide, and indeed the lines where this is the case may correspond to transition lines in 
the phase diagram, since they correspond to lines where the correlation length diverges. It 
is important to note, however, that the high-density isotropic phase is a critical phase, in 
which £ — > oo everywhere in the infinite lattice system, and so the condition that the two 
eigenvalues become degenerate is not a foolproof argument, and must be used with care. A 
standard method for finding phase transition lines in a transfer matrix calculation is to use 
finite size scaling in the form of Nightingale's renormalisation group method[l^|, which is 
based on the scale invariance expected close to critical points for large enough lattice sizes. 
It shows that when solutions exist for the finite-size renormalisation equation 



L L" 



(15) 



then these lines are candidate critical transition lines, although, again, such solutions may 
exist in the high density critical phase without corresponding to transition lines. 

The use of transfer matrices in the determination of the phase diagram is convenient, 
since the partition functions are calculated exactly for infinite strips for any value of the 
parameters given. Since the partition function is known exactly, there are no convergence 
problems, and the full phase diagram may be mapped with relatively little effort. The main 
problem is the size of the matrices, which grow exponentially with the lattice width. This 
strongly limits the maximal width which may be used, here to L = 9. Added to the fact 
that the model contains strong odd/even parity effects, the number of sizes available to more 
advanced finite-size scaling methods is too small to be of much use. In order to be able to use 
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finite-size scaling to calculate critical exponents for various transition lines, we decided to use 
a recently introduced implementation of the corner-transfer-matrix renormalisation group 
(CTMRG) method appropriate for lattice walk models j^]. This method, related to the better 
known density-matrix renormalisation group (DMRG) method, enables the calculation of 
thermodynamic quantities for large lattice sizes, in particular the density of monomers, for 
lattices with the restriction that L x = L y . The large lattice sizes are achieved by iteration 
from smaller lattice sizes. At each iteration the phase space is optimally pruned such that 
the calculation remains within the constraints of the available computer resources and the 
error on the quantities of interest is minimised. For further details on the implementation 



of this model, please see reference 10| and references therein. 



III. RESULTS 

In this section we present the results obtained from the transfer matrix calculations, 
supplemented when necessary with results from the CTMRG method. We start by studying 
the small a regime, where the model is expected to behave like the pure Hydrogen model (a = 
0), studied by Foster and Seno ll| using transfer matrices, and in the Bethe approximation 
by Busano and ftetafl . At the other end of the scale, when a = 1 we recover the pure 
O-point model. Considering values of a close to a = 1 we find new behaviour, not present 
in either the Hydrogen model or the 0-point model. The results found permit the phase 
diagram to be plotted, and finally we will present results for this phase diagram. 

A. Results for a = 0.2 

Foster and Seno 



ll| studied the model for a = using transfer matrices. It was shown 
that the behaviour of this model was different from the standard G-point model. The collapse 
transition was now first order, corresponding to a jump in the density as the low- A' transition 
line is followed. The model also presents a high-A transition line separating the usual 
isotropic dense (liquid) phase from an anisotropic (crystalline) phase. They gave evidence 
that this transition was critical, in contradiction to extended mean-field type calculations 
performed by Buzano and Pretti[3] on the same model, which predicts a first order transition. 
In this section we show that for a = 0.2 we recover a similar behaviour, showing that 
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the crystalline phase and its associated phase transitions persist over a range of values of 
a. Transfer matrix calculations are combined with those of the CTMRG method to obtain 
more accurate results for the critical behaviour of the liquid-crystalline transition, which we 
confirm to be of second order. 

The phase diagram for a = 0.2 is shown in figure HI calculated using transfer matrices. 
In figure H] the \ow-K transition line between the finite-length p = phase and the dense 
phases is found by setting the largest eigenvalue of the transfer matrix to one. When even 
lattice sizes are considered, the largest eigenvalue corresponds to the largest eigenvalue from 
the odd sector of the transfer matrix (A ) for small (3 and the largest eigenvalue of the even 
sector (A e ) when (3 is large. The point where the two cross, A e = A D , is identified as the 
crystallisation transition. This line extends both into the high and low density phases. In 
the low density phase the line corresponds to the crossing of sub-dominant eigenvalues (since 
the largest eigenvalue, corresponding to the empty lattice, is A = 1). This may be identified 
with a disorder line, indicating a change of local order. In the dense phase, however, the 
crossing corresponds to a crossing of the two largest eigenvalues, which from equation (114p 
may be seen to correspond to a divergent correlation length. This line may then be identified 
with a special critical line where the long range order changes. The identification of this line 
with the phase boundary between the liquid and crystalline phases is not straightforward 
since the liquid phase is a critical phase, so it is possible that the line is buried within this 
phase. To verify that this is indeed the transition line, we compare the results with results 
calculated from Nightingale Renormalisation Group arguments (the points shown also in 
figure Hj). In figure M we give finite-size estimates for the transition line for a = 0.2, K = 2, 
showing the coherence of the different methods, and verifying that the eigenvalues cross at 
the transition. 

Transfer matrices are limited by the maximal lattice width that may be obtained, which in 
turn limits the number of lattice sizes which may be used for studying finite-size behaviour. 
For this reason we turn to the CTMRG method, which produces results for large lattice 
sizes, permitting better estimates of the critical exponent, critical density and temperature. 
Unlike the transfer matrix method, where we dealt with infinite strips of finite width, in 
what follows we will be looking at a lattice finite in both directions, with L x = L y = L. 

In order to use finite-size scaling, we consider the scaling form for the density. Here it 
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FIG. 4: Phase diagram for a = 0.2 calculated using eigenvalue crossings (lines), and the phe- 
nomenological RG method (points) for the high-p transition. 
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FIG. 5: Density calculated for a = 0.2, K = 2 using CTMRG with L = 1000 
was convenient to fix K and vary (3, for which we expect the following scaling form: 

Pl{0) = Poo(P) + L^- 2 p(\(3 - (3 c \L l ' v ). (16) 

The scaling form given in equation [16] gives the leading behaviour close to the critical 
point. The corrections to scaling may be expected to be negligible "close enough" to the 
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critical point. It is clear that if the following variables are plotted : 

x = IP-PcIL 1 /", (17) 
y = ( Pi (/3)- Poo )L 2 -^, (18) 

then, close enough to the critical point, the points plotted should fall onto the universal 
curve y = p(x). This phenomenon is known as data collapse. The values of /3 C , v and 
are not known, but can be determined by choosing values which give the best data collapse. 
It is of course useful to know some estimate of (3 C and by some other method, estimates 
which we are able to improve by optimising the collapse of data close to the transition. In 
figure [5] we plot the density as a function of (3 for K = 2 and a lattice L = 1000. It is readily 
seen that p^ « 0.99 and j3 c ~ 0.62 — 0.63. Starting with these initial values we determined 
values of the parameters which gave the best data collapse, and we find: p^ = 0.989 ±0.001, 
/3 C = 0.6222 ± 0.0005 and v = 0.87 ±0.02. The resulting curve y = p(x) is shown in figure [61 
The error bars correspond to the range of values over which the parameters may be varied 
before we clearly no longer have collapse of the data. We limited the lattice sizes to L < 160 
in the study, since, due to the factor L 1 ^ in the variable x, the points which appear in the 
figure are closer to the critical temperature as the lattice size increases. The error in the 
determination of the point is also amplified by the factor L 2 ~ l / V in y. These considerations 
limit the maximum size considered. 

Clearly, in the crystalline phase, the walk will wish to align with one of the lattice 
directions, with a tendency to eject corners from the bulk. Whilst we expect the density of 
corners to differ in the two phases, we note that the data collapse indicates that the density 
of the walk has not reached one, and so the density of corners must be non-zero. This means 
that whilst the fluctuations in the numbers of corners should diverge at the transition, the 
corner density is not a good order parameter. On the other hand the crystalline phase is 
anisotropic, whilst the two other phases are isotropic. A natural order parameter is then 
5p = \p v — ph\, the difference between the densities of vertical and horizontal bonds. In 
the isotropic phases this will vanish, but not in the anisotropic phase. We are not able to 
calculate this quantity with our CTMRG calculation, since the symmetries of the lattice are 
explicitly used in the method^, but we have direct access to this parameter using transfer 
matrices, through 

5p=W\N v (C)-N h (C)\p(C). (19) 
L c 
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FIG. 6: Data collapse for the density close to the transition using data from the CTMRG method 
with a = 0.2, K = 2. The finite-size scaling form of the density is taken with p^ = 0.989, 
j3 c = 0.6222 and v = 0.87 
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FIG. 7: Plot of 5p = \p h - p v \ for a = 0.2 and K = 2. 

This is shown in figure [7| for a = 0.2, K = 2. The peaks of the fluctuations in 5p and the 
density of corners, p c , may be used as estimators for the liquid-crystal phase transition line. 
These are shown, along with other estimates, in figure El 
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FIG. 8: Estimates for (3 C using various methods for a = 0.2 and K = 2. Arrow shows the value 
found using the data collapse method, (3 C = 0.6222. • represent the peaks of the fluctuations of 5p, 
whilst the triangles correspond to the position of the peaks of the fluctuations in the corner density, 
calculated with CTMRG. * gives the position of the solutions to the Nightingale RG method, and 
the squares the position estimated using the condition that A e = A D . 

B. Results for a = 0.8 



In this section we choose to study the phase diagram for a = 0.8, where the model is found 
to have a very different behaviour. The phase diagram calculated using the phenomenological 
renormalisation group is shown in figure EJ Whilst we still find the three phases: the low-i^ 
zero density phase, and the liquid and crystalline phases at higher K, the diagram is quite 
different in appearance. 

As the low-.fr transition line is followed, we have first a type transition, followed later 
by a crystallisation transition. The transition from finite walk to the dense phases is a second 
order transition in the self- avoiding walk class for (3 < (3® becoming first order for f3 > (3®. 
In figure [TP] the corner density fluctuations are plotted along the low-i^ transition line. The 
corner susceptibility is calculated by introducing an additional parameter corresponding to 
a bending energy, and then applying the fluctuation dissipation theorem. Denning /3 CO rn = 
£ C orn/ kT, we define 

Xc = 7^5 • (2°) 



d(3 c 
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FIG. 9: Phase diagram for a = 0.8 calculated using Nightingale's Phenomenological RG method. 
The special transitions along the low-iT line are shown. The first is the 9 point for a = 0.8 whilst 
the second is the collapsed-crystalline transition. The solution of the condition A e = A Q is shown 
for various sizes using points, and can be seen to give a distinct line, which does not converge to 
the collapsed/crystalline transition line. 
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FIG. 10: Fluctuations of the number of corners in the walk for a = 0.8 calculated along the lower 
critical line showing clearly the existence of two transitions. Both transitions appear to be critical. 
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FIG. 11: Bond density calculated for K = 2, a = 0.8 using CTMRG with L = 1000. The circle 
indicates the location of the collapsed/crystalline phase transition calculated using data collapse 
(see figure [T2]) . The transition can be seen to occur at a density p = 1. 

The formation of two peaks may clearly be seen, corresponding to the two special points 
along the low-K line. Interestingly, unlike the a = 0.2 case, the two transitions appear to 
be critical. 

We now investigate the nature of the transition line between the two dense phases. What 
is interesting is that the condition A e = A for even transfer matrices, which coincided 
with this line for a = 0.2, is well within the isotropic collapsed phase here. This is not 
contradictory, since the collapsed phase is critical, and so the condition A e = A Q must 
correspond to a change of order within the critical phase, but not to the phase boundary. 
This is the first indication that the transition line here is different from the transition 
in the previous section. Here again we turn to CTMRG and look for the conditions for 
data collapse. The density, p, appears to have saturated to one (see figure [11]), making 
it impractical to use, however it is expected that the density of corners, p c , should scale 
in the same way, and this is what we use here. The best fit was obtained for p Cj00 = 
0.2744 ± 0.0005, f3 c = 2.349 ± 0.003, and v = 0.96 ± 0.02, and is shown in figure GGl 

In figure [TBI we show different estimates for the critical point for K = 2 and a = 0.8. 
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FIG. 12: Data Collapse of the corner density, p c , for a = 0.8, K = 2 fitted with p c ,oo = 0.2744, (3 C 
2.349, and v = 0.96. 
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FIG. 13: Estimates for (3 C using various methods for a = 0.8, K = 2. Arrow shows the value, 
Pc = 2.349, found using the data collapse method. 

C. The a — f3 phase diagram 

The full phase diagram is expressed in three variables, K, a and /3, and is difficult to 
picture. In this section we present the phase diagram in the a — [3 plane calculated on 
the surface K = K*, where K* is the value of K required for the average length of the 
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FIG. 14: Phase diagram in the a — (3 plane calculated using phenomenological RG and eigenvalue 
crossings. The solid lines show the solutions to the condition A e = A D , expected to coincide with 
the self-avoiding walk-crystalline transition. 

walk to just diverge. This is what is generally calculated in Monte-Carlo simulations, and 
corresponds to the "long polymer in dilute solution" limit. 

When a = 1 the model corresponds to the pure 0-point model, with only one transition 
point on the K = K*(/3) line; the tricritical point. This point is easily found by looking 
at the crossings of the finite-size estimates of v calculated along the line of solutions to 
equation [TH For small (3 these estimates tend to vsaw = 3/4 whilst for large (3 they 
tend to v — 1/2, characteristic of the first order collapsed-walk line in two dimensions. In 
between these two behaviours we find a point, which tends to v® = 4/7. By the way the 
estimates tend to their limiting values, this intermediate point shows up as a crossing in 
the different finite-size estimates. Looking at these estimates as a function of a leads to 
the extended line of tricritical points in the G-point universality class. This is the usual 
method for determining the location of the tricritical point, but it requires the use of three 
lattice widths to determine one estimate. Here we propose a different method. The low-i^ 
transition line is determined by looking for solutions of the phenomenological RG equation 
(|T5|) with the correlation length denned by equation (Ti"4"|) taking Ao = 1 and Ai is the largest 
eigenvalue taken from the odd or even sectors of the transfer matrix. A tricritical point has an 
additional correlation length which diverges, corresponding to the two relevant directions in 
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FIG. 15: Phase diagram in the a — j3 plane calculated using the peaks of xp an d X&p calculated 
using transfer matrices. The peaks of p fluctuations pick out the line of G points, whilst the peak 
of bp fluctuations pick out the transition between the isotropic collapsed phase and the anisotropic 
crystalline phase. These two lines merge to form the first order transition line separating the SAW 
phase from the crystalline phase. 



the renormalisation group sense. We look for the solutions of equation (1151) with a correlation 
length calculated using the largest eigenvalue from the odd and even sectors. This method 
only requires two lattice widths to estimate the location of the tricritical point. Additionally, 
the high-.fr transition line is found using (jl~5l) with these same two eigenvalues. The method 
described therefore also locates the position of the crystallisation transition along the low-i^ 
transition line. The phase diagram calculated by this method is shown in figure [HI For 
small a there is no solution. However this region of the phase diagram corresponds to the 
region which is expected to behave like the pure hydrogen-bonding model, and therefore the 
condition A e = A Q corresponds to the first order transition line. 

In figure [15] the phase diagram is calculated using the peaks of the two susceptibilities. 
The upper line is calculated by looking for the peak of the fluctuations of 5p whilst the lower 
line is calculated using the peaks of the fluctuations of p. For a less than about 0.4 the two 
sets of lines merge and give estimates for the single hydrogen-bonding like first order line, 
whilst for values of a larger than about 0.5 the two sets of lines are distinct, the lower line 
corresponding to the line of like tricritical points, whilst the upper line corresponds to the 
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critical crystallisation line. Somewhere in the region a = 0.4 — > 0.5 these two lines merge 
into a higher order multi-critical point. 

Looking closely at figure [T41 it may be seen that the upper and lower transition lines tend 
to come together in the region a = 0.3 — ► 0.5, with the lower line developing a plateau. 
It is probable that the multicritical point is not located at the cusp where solutions end, 
but at a higher value of a, probably in the same range of values. The methods employed 
in this article were not able to determine this point more accurately. The CTMRG, which 
enables larger sizes to be obtained, becomes impractical in this region, particularly when 
the number of constraints required to define the point is considered. 



IV. DISCUSSION 



There are now a number of similar models which display an anisotropic crystalline phase, 
with a variety of different types of high-density transition. The first is the vertex-interacting 
self-avoiding walk due to Blote and Nienhuisj^j], which displays an Ising like high-density 
transition with u = 1. The Hydrogen-bonding model was shown to also have a critical 



transition 



11], but in a different universality class, which is confirmed in this article, where 



the value v ~ 0.87 is found. Lastly the bond-interacting O-point model|8i] which is conjec- 
tured to have a softer, higher order critical transition. 

What is interesting in the model presented here is that two different high-.fr critical 
behaviours are displayed in one model. For smaller values of a the transition from the 
collapsed to crystalline phase is of the Hydrogen- Bonding class, and occurs at densities 
which are close to p = 1 but on close inspection we clearly have p < 1, as may be seen in 
figure [51 

Data collapse for a = 0.8 gave v = 0.96 for the best fit, however the fitting was less clear 
than for a = 0.2 and it is possible that the correction terms are more important. The value 
of v calculated would then be an effective exponent. It is tempting to conjecture that the 
true value of v = 1, in analogy with the vertex-interacting model. However, as may be seen 
in figure [11] the transition may be seen to occur well after the density saturates to p = 1. If 
this is the case, then the walk is essentially a Hamiltonian walk at the transition, and looks 
very much like the model with a penalty for the formation of corners in the Hamiltonian walk 
limit studied by Saleurjl5], where he conjectured that the transition was of infinite order 
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BKT transition of the same type as in the F-model, but is in contradiction to the results 
presented here. This contradiction was first seen in the Hydrogen model in the Hamiltonian 
limit 11], where the transfer matrix results also gave estimates close to v = 1. This is a 



point which warrants further investigation. 

The results found here are to some extent confirmed by a Monte-Carlo study, mainly 
in three dimensions, which has appeared during the final stages of this work[3]. Notably, 
based on a flat-PERM study, a similar phase transition to that presented in figures [bil and [T5 
is found. The transition from the collapsed phase to the crystalline phase along the surface 
where the walk length just diverges was seen to be probably critical, as is the case here. 
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